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ABSTRACT 

The rapidly oscillating Ap stars are of importance for studying the atmospheric structure of 
stars where the process of chemical element diffusion is significant. We have performed a 
survey for rapid oscillations in a sample of 9 luminous Ap stars, selected from their location 
in the colour-magnitude diagram as more evolved main-sequence Ap stars that are inside the 
instability strip for rapidly oscillating Ap (roAp) stars. Until recently this region was devoid 
of stars with observed rapid pulsations. We used the VLT UV- Visual Echelle Spectrograph 
(UVES) to obtain high time resolution spectroscopy to make the first systematic spectro- 
scopic search for rapid oscillations in this region of the roAp instability strip. We report 9 
null-detections with upper limits for radial-velocity amplitudes of 20 — 65 ms -1 and preci- 
sions of a = 7 — 20 ms _1 for combinations of Nd and Pr lines. Cross-correlations confirm 
these null-results. At least six stars are magnetic and we provide magnetic field measurements 
for four of them, of which three are newly discovered magnetic stars. It is found that four stars 
have magnetic fields smaller than ~ 2 kG, which according to theoretical predictions might be 
insufficient for suppressing envelope convection around the magnetic poles for more evolved 
Ap stars. Suppression of convection is expected to be essential for the opacity mechanism act- 
ing in the hydrogen ionisation zone to drive the high-overtone roAp pulsations efficiently. Our 
null-results suggest that the more evolved roAp stars may require particularly strong magnetic 
fields to pulsate. Three of the studied stars do, however, have magnetic fields stronger than 
5kG. 

Key words: stars: individual: HD 107107; HD 110072; HD 131750; HD 132322; 
HD 151301; HD 170565; HD 197417; HD 204367; HD 208217- stars: magnetic fields -stars: 
oscillations - stars: variables: other 



1 INTRODUCTION 

Why do some stars oscillate while others do not? This question is 
particularly relevant for the class of rapidly oscillating Ap (roAp) 
stars for w hich the disc overy of high-frequency oscillations came as 
a surprise (Kurtz 1982). The <5Scuti pulsations of stars in this area 
of the colour-magnitude diagram (CMD) were not theoretically 
predicted to produce the observed high-frequency, high-overtone 
modes. Subsequently, considerable observational and theoretical 
efforts have been able to describe most pulsati on propertie s of 
roAp stars through i) the oblique pulsator model (Kurtz 1 9821). ii) 



the pr obable driving mechanism of the pulsations (Balmforth et al. 
l200lh and iii) their link with the m agnetic fields present in Ap stars 
(see, e.g.. lCunhall200rj : ISaioll2005l) . 



* Based on observations collected at the European Southern Observatory, 
Paranal, Chile, as part of programmes 075.D-0145 (A), 078.D-0080(A), 
072.D-0138(A) and 077.D-0150(A). 
f E-mail: lmfreyhammer @ uclan.ac.uk 



Ap stars are chemically peculiar stars that range from early- 
B to early-F spectral types with the majority having detectable 
magnetic fields with strengths of a few 10 2 to a few 10 4 G 
(Bvchkov, Bychkova, & Madei 20031, according to whom ~ 55 per 
cent have mean longitudinal fields above 400 G). For the cool Ap 
stars (also called CP2 stars), abundance anomalies are typically 
visible in their spectra as abnormally strong absorption lines, par- 
ticularly of rare-earth elements (REEs). Candidate roAp stars can 
therefore be selected photometrically due to the influence of these 
spectral features on narrow-ba nd photometric i ndices such as those 
of the Stromgren filter system {Martinez 1993, see also Sect. |2.1| >. 

The high-overtone roAp oscillations are thought to be driven 
by th e opacity-mechanism acting in the partial hydrogen ionisation 
zone toziembowski & Gooddl 19961 ; iBalmforth et allhoOll ; ICunhal 
120021 : Saio 2005). The frequency range of the excited modes is re- 
lated to the presence of strong magnetic fields in roAp stars that i) 
directly stabilise the low-frequency, low-order S Scuti oscillations, 
excite d by the ^ -mechanism acting in the Hell partial ionisation 
zone Jsaioll2005h , and ii) indirectly enhance the driving of high- 
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order oscillations by t he w-mechanism acti ng in the partial hydro- 
gen ionisation zone tealmforth e"t~al]|200lh . In particular, the sta- 
bilisation of the 8 Scut i modes ca n be explained by the dissipation 
of slow Alfven waves (Saio 2005) and by the suppression of turbu- 
lent motion in the outer layers of the magnetic pole regions of roAp 
stars, which speeds up gravitat ional settling of hel ium and drains it 
from the He II ionisation zone dTheado et al.l2005T) . An implication 
of the high radial order of the observed modes of roAp stars is that 
the observable pulsation amplitudes become depth dependent. This 
pulsational structure provides a unique opportunity for studying the 
magneto-acoustic structure of the pulsations in 3D because of the 
element stratification in the atmospheres of these stars. Different el- 
ements, such as Fe, Pr, Nd, Eu and H, act as tracers for the pulsation 
structure due to their different properties regarding formation depth 
and extent of formation regions. Examples of dept h-dependent 
measu r ements of pulsations are now numerous, such as Kurtz et al. 
j2005b : lEikin. Kurtz. & Mafhvsl d2005h ; lRvabchikova et alj IXXm . 
The roAp stars are thus unique objects for studying the interactions 
among stellar pulsation, magnetic fields and atomic diffusion. The 
latter is important for, e.g., estimates of globular cluster ages where 
He settling is significant, pulsation driving mechanisms in sdB and 
[3 Cephei stars where radiative levitation of Fe is needed for insta- 
bility, and the Standard Solar Model where both He settling and 
radiative levitation of some metals are included. 

Only 37 roAp stars are know n at present dKurtz et"al] [2006b ; 
iTiwari. Chaubev, & Pandevl 120071) despite seve ral searches for 
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mostly based on photometric indices pointing to chemical peculiar- 
ities in the spectra. There is at present no clear correlation between 
radial velocity (RV) a mplitude and phot ometric amplitude in roAp 
stars (see Table 1 of iKurtz et alj|2006bl) . Fo r a (non-exhaustive ) 
list of spectro scopic studies of roAp stars, see lKurtz et alj J2006al) . 
Cunha ( 2002) calculated a theoretical instability strip for roAp stars 
and compared it to locations of 16 known roAp stars. The region 
around the terminal age main-sequence is remarkably devoid of 
pulsators which, as Cunha points out, could be an observational 
bias. She predicted that more luminous and evolved Ap stars may 
pulsate with lower frequencies (0.67 — 0.83 mHz) which makes it 
harder to detect them in the typically short high-speed observing 
runs. Note that this frequency range overlaps with the highest 8 Se t 
frequencies, such as for HD 34282 (0.92 mHz. lAmado et al]2004l) . 
Alternatively, she proposed that this absence of known evolved 
roAp stars could reflect a real deficiency in the fraction of pul- 
sators in that region of the HR diagram, resulting from the fact that 
the magnetic field is less likely to suppress envel ope convection in 
more evolved stars. The photometric survey by Martinez & Kurtz 
dl994l) was indeed less sensitive to low-amplitude roAp periods 
longer than 15 min and got a null-result for HD 1 161 14 for which 
a 21-min oscillat ion was recently detected spectroscopically by 
Elki n et alj d2005t) . A dedicated survey of luminous Ap stars there- 
fore seems pertinent. 

Ap stars inside the roAp instability strip that do not exhibit 
any detectable variability are called non-oscillatin g Ap (noAp) 
stars. A s seen in, e.g., the astrometric HR-diagram by Hubr ig et al.1 
d2005l . their figure 2) for roAp and Ap stars, the apparent noAp 
stars occupy essentially the same regions as the roAp stars. How- 
ever, the noAp star s appear to be systematically more evolved 
than the roAp star s dNorth et al.lll997l ; lHandler & Paunzen|[l999l ; 
iHubrig et alj|2000l) . To fully understand the mechanism responsi- 



ble for driving oscillations in roAp stars, it is essential to confirm 
and understand why so many stars with similar characteristics are 
apparently stable against pulsations. The apparent absence of os- 
cillations in noAp stars does not, however, necessarily mean that 
oscillations are suppressed. For instance, detection of oscillations 
may depen d on lifetimes of excited modes (may be days for some 
roAp stars, Handler 2004), on beating between multiple modes, on 
the geometry of the magnetic field's orientation and the aspect in 
which the stellar surface is observed together with the surface dis- 
tribution of the chemical elements used to detect the oscillations. 
Further, the signal-to-noise ratio (S/N) of the obtained photome- 
try or spectroscopy may simply be too low to detect the typically 
low-amplitude roAp pulsations. Finding high precision null-results 
for a statistically significant sample of roAp candidates would be a 
safe basis for concluding whether noAp stars exist or whether the 
current lack of known luminous roAp stars is an observational bias. 
From a theoretical point of view, stars located inside the theoreti- 
cal instability strip may be noAp stars if the magnetic field is too 
weak to effectively suppress convection in the stellar envelope. A 
good test sample should, therefore, also include stars with strong 
magnetic fields. 

To try to answer the question of whether the absence of ob- 
served variability among the luminous Ap stars is an observational 
selection effect only, or in fact due to intrinsic properties, we have 
studied 9 such stars at high radial-velocity precision. Our immedi- 
ate goal with UVES was to search these luminous Ap stars for roAp 
oscillations at high radial velocity precision. Table[T]lists the known 
properties of the 9 targets and gives an observing log for the col- 
lected spectra. We find no pulsation in any of the 9 stars; Tables|4] 
and[5]below present the null-results of the frequency analyses. 

In the following sections, we describe selection and observa- 
tion of the targets along with the data reduction in Sect. 2, then 
follows the data analyses including estimation of physical param- 
eters and the radial-velocity analysis in Sect. 3. We show that our 
analyses reach the same precision as other radial-velocity studies 
of roAp stars in the literature and we then give the results for each 
target on a star-by-star basis. Finally we discuss the null-results in 
Sect. 4. 



2 SELECTION OF THE CANDIDATES, OBSERVATIONS 
AND DATA REDUCTION 

2.1 Selection 

We selected 9 stars based on Stromg ren and /3 indices from the 
Cape catalogue of cool Ap stars bv lMartinezl dl993l) . and luminosi- 
ties from Hipparcos parallaxes. For Stromgren photometry, the ci 
index is normally an indicator of luminosity while mi is an indi- 
cator of metallicity. The /3 index is sensitive to temperatures for 
stars in the range A3-F2. Also b — y indicates temperature, but is 
influenced by reddening. The relative indices Sci and 5mi indi- 
cate how a given star deviates i n ci and mi from 'normal' zero-age 
main sequence stars (see also ICrawford|[l979l) . A more negative 
value of <5mi indicates stronger metallicity. For normal stars <5ci is 
a luminosity indicator, while for Ap stars it is depressed by heavy 
line blocking and hence is not a reliable indicator of luminosity. A 
negative 8c\ is, in fact, indicative of strong peculiarity. However, 
because ci increases with luminosity, luminous Ap stars may have 
positive Sci which, without independent information on the lumi- 
nosity, makes their chemical composition appear 'normal' in this 
index. A negative 5c\ has frequently been used for target selection 
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Table 1. Observing log indicating number of observations and the mid-series heliocentric Julian Date for each series of 
spectra. Exposure times were 40 s, except for the fainter stars HD 1 10072 and HD 170565 where 80 s was used; readout and 
setup times were 24 — 27 s. The last two columns indicate S/N measured in the continuum based on the standard deviation 
in the residual flux from two consecutive spectra. Lower and Upper respectively refer to wavelengths below or above the gap 
between the two CCDs at 6000 A, while the ranges in parentheses correspond to the S/N range for all spectra in that series. 
The S/N estimates assume random noise only. 



Star 


#spec 


Q2000.0 


feooo.o 


HJDobs 


Time 


S/N Lower 


S/N Upper 






(h:m:s) 


(°:V) 


(d) 


(h) 






HD 107107 


111 


12:19:04.8 


-40:09:47 


2453509.504 


1.94 


79 (77-82) 


62 (56-77) 


HD 1 10072 


69 


12:39:50.2 


-34:22:30 


2453509.591 


2.06 


50 (43-55) 


38 (32-49) 


HD 131750 


111 


14:56:20.7 


-30:52:36 


2453509.681 


1.98 


80 (74-87) 


66 (57-74) 


HD 132322 


111 


15:01:36.0 


-63:55:35 


2453510.519 


2.03 


151 (147-155) 


141 (132-152) 


HD 151301 


111 


16:49:28.3 


-54:26:47 


2453510.715 


1.95 


82 (76-84) 


67 (60-79) 


HD 170565 


85 


18:30:08.2 


-02:35:26 


2453509.778 


2.50 


89 (83-94) 


77 (69-86) 


HD 197417 


125 


20:48:48.1 


-72:12:44 


2453509.879 


2.20 


107 (94-113) 


93 (79-110) 


HD 204367 


111 


21:28:41.2 


-25:38:39 


2453510.798 


1.96 


113 (108-120) 


101 (93-115) 


HD 208217 


138 


21:56:56.4 


-61:50:44 


2453510.896 


2.52 


156 (142-165) 


148 (128-159) 



in previous searches for roAp stars; as a consequence, luminous Ap 
stars have seldom been tested for pulsation. Our sample of 9 roAp 
candidates was therefore hand-picked among Ap stars with Hippar- 
cos parallaxes; their main physical properties are given in Table[2] 
Fig.Q]shows that the astrometric luminosities (Table[2] column 11) 
place the chosen sample among the more evolved Ap stars (towards 
the end of their core-hydrogen burning phase). The corresponding 
temperatures (column 9) are described in Sect. |3.1.Tl The absolute 
magnitude in the V-band and the stellar luminosity were derived 
using the standard relations: 



M v = m v + 5 + 5 log ty - A x 



(1) 



where the trigonometric parallax ir is measured in arc seconds and 
the interstellar extinction in the V-band is Ay = 4.3i5(& — y) = 
4.3 • 0.74 E{B - V), and 



L Mv+BC- M boli0 
log — = 



L 







2.5 



(2) 



where we used M bol q = 4.72. Column 8 in T able[2]gives redden- 
ing read from the maps bv lBurstein & Heilesl dl982h and columns 
11 and 12 give the luminosities derived with and without redden- 
ing. The reddening maps have errors in E(B — V) of 0.01 mag 
or 10 per cent (whichever is the greatest) and resolution of 0.03 
mag. HD 170565, HD 132322 and HD 151301 are near the galactic 
plane and outside these maps so a lower limit on reddening was 
used. These reddening estimates are in all cases conservative val- 
ues but contribute little to the luminosities where errors are dom- 
inated by the errors from the parallaxes. The only exception is 
HD 132322 which has a p recise distance. No Lutz-Kelker correc- 
tion dLutz & KelkeJ 19730 was applied to the luminosities. The val- 
ues for projected rotational velocities and magnetic fields in Table|2] 
are described further in Sects. [3~TTT| and [3~1.2| 

Fig.Q] shows the locations of these stars in a colour- 
m agnitude diagram (CMP) sup erposed with evolutionary tracks 
by IChristensen- Dalsgaard il993h and locations indicated for stars 
that a re predicted to be pulsationally stable or unstable dCunhal 
l2002h . Lum inosities are from distances based on Hipparcos mea- 
surements (ESa| 1 19971) with considerable errors from the par- 
all axes indicated. Temperatur es were estimated from the grids 
bv | Moon & Dworetskvl d 19851) based on [3 photometry dMartinezI 
ll993l) alone.~ 



2.2 Observations and data reduction 

Our spectroscopic observations were obtained at the Very Large 
Telescope (VLT) using the Ultraviolet and Visual Echelle Spectro- 
graph (UVES). The fast readout mode of 625 kpixs" 1 was used 
for reading out the two CCD-mosaics on the red arm of UVES in 
the 4-port low gain mode. We observed our 9 targets about 2 h each 
(see the observing log in Table 1) on the two nights of 18 — 19 
May and 19 - 20 May 2005 (JD2453509 and JD2453510), using 
exposure times of 40 s with 24 — 27 s readout and overhead times, 
giving a mean time resolution of 64 s. However, for the faintest 
stars HD 110072 and HD 170565, exposure times were doubled to 
80 s, which reduced the time resolution to 107 s. The total number 
of spectra per target ranges from 69 to 138. To ensure high count 
rates, an image slicer was used. The spectra were processed with 
the UVES pipeline and ESO MIDAS package to extract ID spectra 
using nightly Hatfield spectra and thorium-argon calibrations ob- 
tained at each telescope pointing. The spectra were extracted us- 
ing the 'average extraction' option and the flatfielded spectra from 
the two CCDs on the mosaic were rebinned to same step size in 
wavelength and combined in one spectrum. The rectification was 
performed in 3 steps. First, all spectra from the same night were 
normalised with a spline-fit to a high-S/N spectrum of that night, 
only considering large-scale patterns, such as slopes. Next an un- 
dulating continuum pattern that followed the spectral orders was 
fitted with splines using averaged spectra of all spectra each night 
and then applied to eliminate the undulations in individual spectra. 
Finally, all spectra of the same series (i.e. of the same star) were rec- 
tified relative to a mean spectrum of that series. Removal of effects 
from cosmic rays and CCD blemishes in the spectra was made in 
two steps: cosmic rays were identified and removed by comparing 
each 2D CCD image with the previous and subsequent one, and the 
final, rectified ID spectrum was corrected by a routine that identi- 
fied and removed sharp emission features of non-stellar origin. For 
all spectra, the region AA 6515 — 6535 A suffers from spectrum-to- 
spectrum variability on the level of 10 per cent. A gap occurs in the 
region AA 5963 — 6032 A and is caused by the gap between the two 
CCD mosaic halves (referred to as 'lower' and 'upper'). 

For unknown reasons, our UVES reductions resulted in an in- 
creased noise in the spectra redwards of the gap, which varies from 
star to star but is particularly pronounced for HD 1 10072. This af- 
fects the accuracy of v sin i, magnetic and radial-velocity measure- 
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Table 2. Main properties of the 9 luminous Ap stars. The columns give: ( 2 — 6) Stromgren indice s from Martinez 1 1993), (7) the parallax it from Hipparcos - 
the values in columns 10 - 12 are based on this value, (8) reddening from lBurstein & HeilesH 19821) fixed to 0.21 for HD 170565, HD 132322 and HD 151301, 
(9) temperatures derived from the (3 index, (10, 12) absolute magnitu de and lumino sity assuming no reddening, while luminosity in Col. 1 1 accounts for 
reddening. Bolometric correction is from interpolation in the tables of Flower 1 1996). Projected rotation velocities vsmi (13) with estimated errors were 
obtained with the measurements of the mean quadratic magnetic field (14). Magnetic field measurements from the literature are specified with superscripts. 
The numbers in the last column (15) associate the stars with labels in Fig.[T] 



Star 


V 


b-y 




P 


7T 


E B _y 




My 


log L/Lq 


log L/L 


vsini 


(B q ) 


Note 


HD 


mag 


mag 


mag mag 


mag 


(mas) 


mag 


K 


mag 






kms — 1 


kG 




107107 


8.734 


0.074 


0.234 0.835 


2.864 3.42 ± 1.11 


0.10-0.14 


8300 


1.404 


1.33 


1.68 


10.5 ±0.5 


5.2 ±0.4 


#5 


110072 


10.104 


0.227 


0.281 0.487 


2.759 


2.45 ± 1.44 


0.06-0.09 7300 


2.050 


1.07 


1.25 


3.3 ± 0.5 


1.5 ±0.6 


#9 


131750 


8.549 


0.121 


0.260 0.710 


2.847 


2.86 ± 1.21 


0.06-0.12 


8100 


0.831 


1.56 


1.79 


25.3 ± 1.8 


5.3 ±3.3 


#2 


132322 


7.357 


0.129 


0.227 0.943 


2.894 


5.78 ± 0.73 


>0.21 


8600 


1.167 


1.44 


2.03 


34.3 ± 1.8 


^ 6.0 a 


#3 


151301 


8.541 


0.222 


0.244 0.666 


2.827 


3.56 ± 1.40 


>0.21 


8000 


1.298 


1.37 


1.90 


13.7 ± 1.1 


^ 2.4 


#7 


170565 


9.130 


0.269 


0.253 0.628 


2.825 


2.72 ± 1.45 


>0.21 


7900 


1.303 


1.37 


1.90 


18.0 ±3.0 


N/A 6 


#6 


197417 


8.009 


0.038 


0.233 0.939 


2.877 


3.23 ±0.84 


0.03-0.06 


8400 


0.555 


1.68 


1.82 


25.5 ± 0.4 


^ 2.2 


#1 


204367 


7.830 


0.039 


0.220 0.996 


2.899 


5.09 ±0.92 


0.03-0.06 


8700 


1.364 


1.36 


1.49 


10.8 ± 0.3 


0.4 


#4 


208217 


7.199 


0.102 


0.274 0.626 


2.816 


6.83 ±0.90 


< 0.03 


7700 


1.371 


1.34 


1.37 


10.8 ± 0.7 


8.0±0.5 C 


#8 



a {B z ) = 0.36 ± 0.05 kG (Hubrig etal. 2006) . b (B r ) = 1.76 ±0.17 kG fKudrvavtsev et al. 2006) . 
C (B) = 7.96 ± 0.59 kG (Mathvs et al. 1997) 




Figure 1. Colour magnitude diagram with the 9 UVES targets superposed and labelled as in Tableg] #1 (HD 197417), #2 (HD 131750), #3 
(HP 132322), #4 (HP 204367 ), #5 (HP 107107), #6 (HP 170565), #7 (HP 151301), #8 (HP 208217) and #9 (HP 110072). Stellar evolution tracks are from 
IChriste nsen-Palsgaard 1 1993) and circles indicate theoretically predicted locations for stars with pulsationally stable (filled) and unstable (open) modes. 
Luminosities are based on Hipparcos parallaxes and the associated uncertainties are indicated with error bars. Predicted pulsation frequencies are given in 
mHz. 
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ments by increasing the noise for this spectral region, in some cases 
by 33 per cent. The effect is most clear when comparing the noise 
of cross-correlation results in the spectral regions bluer and redder 
than the gap (Table[5}, or when comparing noise in co-added series 
of spectra of individual stars on both sides of the gap. The effect 
varies for individual stars and we suspect the origin of the problem 
to be found in the reduction of the spectra. During observations, it 
is more difficult to keep an object on the sheer when observing very 
close to zenith. This also affects the S/N and was indeed the case 
for HD 110072 and HD 131750 that were both observed at about 
10 degree distance from zenith. However, as the performed radial 
velocity analyses utilise regions on both sides of the gap, includ- 
ing the unaffected 1000 A wide region below the gap, and do reach 
precisions comparable to those published in the literature for roAp 
stars, the current quality of the reductions suffice for the aims of 
this investigation. 

The zero point of the absolute wavelength calibration is of lit- 
tle importance for the differential radial-velocity analysis, but was 
found accurate t o the level of 300 ± 200 ms -1 based on the tel- 
luric line list by iGriffin & Griffin] dl973l) . This is as good as one 
may expect from radial velocities of telluric lines that depend on 
conditions in the Earth's upper atmosphere. Barycentric velocity 
corrections were recorded in the FITS headers of the reduced spec- 
tra for the later analysis of the velocity fields, and not included in 
the wavelength calibration. 



Table 3. List of laboratory wavelengths of the most important lines used in 
the radial-velocity measurements based on the Centre-of-Gravity method. 
Note that due to broadening and abundance differences among the roAp 
candidates, not all lines could be used for all stars, and in some cases line 
blending made it necessary to include close lines in the measurements. Lab- 
oratory wavelengths have been taken from NIST, DREAM and VALD. For 
comparison with the stellar lines a few telluric lines were used (A 5919.2, 
5946.0, 5949.2, 6278.9 and 6889.9 A). 



Line Wavelength (A) 

Hi 6562.852 

Li I 6707.761 

Nal 5895.924,5889.951 

Cal 6122.217,6162.173 

Cm 5280.054,5310.687 

Fe I 5434.524, 6136.614, 6137.694 

Fell 5061.718,5284.109 

Ball 5853.675,6141.713 

Cell 5077.854, 5117.946, 5147.565, 5459.193, 5468.371, 

5613.694, 5680.261, 5711.437, 5858.546 
Prll 5220.108,5605.642,6114.381 

Prill 5208.507, 5284.693, 5299.993, 5844.408, 5956.043, 6053.003, 

6090.010, 6160.233, 6195.619, 6866.793 
Ndll 5319.811, 5356.967, 5361.467, 5620.594, 6425.779, 6550.178 
Ndm 5085.001, 5203.902, 5286.764, 5825.857, 5845.068, 

6145.072, 6327.244, 6550.326 
Eull 6437.680,6645.064 
Gdll 5860.727 



3 DATA ANALYSIS AND RESULTS 

Our aim with this study was to search the spectra of nine roAp can- 
didates for rapid oscillations and through simple estimates of phys- 
ical properties to characterise the detected noAp or roAp stars. For 
this purpose, we used the cross-correlation technique and centre-of- 
gravity procedure to search for rapid pulsation using, respectively, 
Doppler shifts of whole wavelength regions and of individual and 
combined line profiles. Stellar properties were estimated with pho- 
tometry and/or by comparison of observed spectra to synthetic line 
profiles. 

For the purpose of spectral line identificat ion, a synth etic com- 
parison spectrum was produced with SYNTH JPiskunovll 19921) us- 
ing a Kurucz stellar atmosphere model for T e s = 8750 K, turbu- 
lent velocity 2kms~ 1 , log g — 4.0 and slightly increased solar 
metal abundance \ogNz/N z q = +0.5. Atomic lin e data were 
taken from the Vienna Atomic Line Database (VALD. lKupka et all 
fl999T) for increased abundances mainly for the elements Nd, Pr, 
Sr, Cr, Eu, as given in Table lAl l Other sources used for line data 
were the atomic database NISTy and the Database on Rare Earth 
Elements at Mons University (DREAV0) through its implementa- 
tion in the VALD. This model's temperature is less optimal for the 
coolest of the studied stars, but still appropriate for selecting lines 
for the radial-velocity analyses. 

In the following section we describe the estimation of stellar 
physical parameters for the examined sample, the radial-velocity 
measurements and frequency analyses, and a performance test car- 
ried out on UVES data for two known roAp stars. 



3.1 Analysis 

The present study is a survey for variability, and therefore an anal- 
ysis of abundances, temperatures, surface gravities and projected 
rotational velocities is outside the scope of this paper. Such an 
analysis must take the chemical element stratification, and possi- 
bly also temperature gradients, into account. We do, however, in 
Table|2]give estimates of vsini, effective temperatures and mag- 
netic fields of the studied sample. 

3.1.1 Temperatures and rotational velocities 

The spectra of Ap stars have flux distributions that are deformed 
by strongly peculiar elemental abundances. Calibrations based on 
Stromgren indices may therefore not provide reliable estimates of 
temperatures and surface gravity, log g. The (3 index, however, re- 
mains largely unaffected by this. Effective temperatures wer e thus 
estimated with the Co, P grids b y [Moon & Dworetskvl |l985) using 
(3 from Tablel2l (from lMartinezll 1993 ). We allowed the value of Co 



1 http://physics.nist.gov 

2 http://w3.umh.ac.be/~astro/dream.shtml 



to be free in these grids, fixed (3 and assumed 3.5 ^ log g ^ 4.0. 
The resulting photometric temperatures are given in Table[2] An 
additional check on T e s was made by comparing synthetic line 
profiles to the observed spectra. This method is sensitive to normal- 
isation errors in the Ha region (the above -mentioned undula tions), 
and the core-wing anomaly in this line dCowlev et al . 2001) that 
prohibits fitting the atomic line profile with normal models. How- 
ever, the average deviation of the photometric and spectroscopic 
temperatures is only 5 per cent, with the exception of HD 131750 
which has a 17 per cent higher spectroscopic temperature. Moon & 
Dworetsky's grids are furthermore ambiguous for T c g > 8500 K 
and in those cases the spectroscopic temperatures of HD 132322 
and HD 204367 were 8250 - 9250 K and 8000 - 8250 K, re- 
spectively. The synthetic line profil es were based on Kurucz at- 
mospheres dCastelli & Kuruczll2003l) and atomic line data from the 
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Figure 2. Comparison of the Ha profiles of all examined Ap stars. To increase the readability, the individual spectra have been offset in flux, and appear with 
increasing temperature upwards. Note the secondary spectrum superposed on that of HD 1 10072, and the clearly double Nd III 6550.32 A lines for HD 170565 
and HD 151301 (indicative of abundance spots). Each spectrum was made by co-adding all spectra in the series for that star. The spectrum of HD 110072 was 
smoothed over two pixels. 



VALD database using increased REE abundances (identical for all 
models). The grid of models had T cff = 7000 - 10500 K in steps of 
500 K with log g — 3.5, 4.0 and, in a few cases, 4.5. Best fits were 
obtained for log g — 3.5 — 4.0, mainly to the lower side of this 
range. Better estimates of T a s require a detailed abundance analy- 
sis that takes stratification and magnetic fields into account, which 
will not be done here. 

Projected rotation velocity estimates were made by measuring 
the Fe I lines A 5434.52 and 5576.08 A that are rather insensitive 
to magnetic broadening. Synthetic models for log g = 3.5 were 
compared to the spectra for a range of iron abundances and rota- 
tional broadening, using models corresponding in temperature to 
the individual stars. The macroturbulence was varied in the range 
1 — 4 kms -1 . When strong line blending hampered the measure- 
ments, additional iron lines were used to constrain the estimates. 
The derived velocities were then refined with synthetic models that 
took magnetic broadening into account (SYNTHMAG) jPiskunovl 
1999). The analyses of several, typically 20-30, iron lines that de- 
termine the m ean quadratic magnetic fi elds also give precise rota- 
tion velocities dMathvs & H ubrig 2006). These values are formally 
more precise than, but consistent with, the above estimates and 
are therefore given in Table[2] except for the case of HD 170565 



where only the first estimate is given. It should be noted, though, 
that the values of the rotation velocities that are obtained as part of 
the analysis performed to determine the mean quadratic magnetic 
field are upper limits on the v sin i, since they may actually in- 
clude contributions from other line broadening effects that are pro- 
portional to wavelength, such as micro- or macroturbulence. (In- 
strumental and thermal broadening are however duly isolated - see 
iMathvs & Hubrid J2006h for details.) However these contributions 
are mostly negligible, except in the slowest-rotating stars. 



3.1.2 Magnetic fields 

Known roAp stars have strong magnetic fields dKurtz et al.l2006ah . 
and we therefore searched for magnetically resolved or broadened 
lines in the spectra of the roAp candidates. Many Ap stars have 
fields that are measurable only with polarimetry. It takes a strong 
field, typically exceeding ~ 1.5 kG, combined with a slow pro- 
jected rotation rate (smaller than v sin i ~ 10 kms~ ) to produce 
magn etically resolved lines by the Zeeman effect (Math vs~et al.l 
1 1997b . 

In the simplest cases of spectral lines corresponding to doublet 
or triplet Zeeman patterns, simple formulae can be applied to deter- 
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Table 4. Frequency analyses from combination of individual line shifts determined for multiple lines of different species with the Centre-of-Gravity method. 
The highest amplitude in each stacked amplitude spectrum is given as A max . 'n' indicates number of measurements, typically number of lines or components 
in case of double structures. Linear shifts have been fitted and removed separately from each time series. See TablefJJfor a listing of the lines typically used. 
Low-frequency peaks (0.4 mHz or below) were excluded from the noise estimates and not considered for significance. 'Nd,Pr' indicates combined radial 
velocity series for all possible Nd and Pr lines. 



Star Nd II Nd III Ce II Pr II Pr III Nd.Pr Eu II Telluric H core 





a A, 


Tiax ft 


a A_ 


tnax H 


a A 


max 


n 


a A 


max ft 


a A, 


max n 


G A ma x n 


a A 


max n 


(7 A ma x n 


a A 


max 


HD 


( m s 


- 1 ) 


( m s 


- 1 ) 


( m s" 






( m s 


- 1 ) 


( m s 


- 1 ) 


( m s' 


- 1 ) 


( m s 


- 1 ) 


( m s" 




( m s — 


x ) 


107107 


28 


82 4 


28 


108 6 


32 


118 


6 


44 


133 2 


16 


59 10 


12 


39 25 


23 


88 2 


20 


62 2 


20 


65 


110072 


27 


92 5 


16 


52 8 


30 


100 


5 


31 


110 5 


15 


55 10 


10 


42 28 


31 


95 2 


20 


57 2 


52 


135 


131750 


38 


132 4 


30 


102 5 


40 


150 


6 


59 


177 2 


29 


112 12 


18 


58 23 


48 


161 2 


20 


57 3 


21 


72 


132322 


32 


90 5 


20 


69 11 


N/A 


N/A 





58 


196 1 


22 


70 8 


13 


37 26 


22 


72 4 


10 


30 3 


19 


59 


151301 


24 


84 3 


11 


35 9 


62 


237 


1 


41 


136 2 


12 


39 9 


8 


24 25 


20 


86 4 


13 


42 4 


19 


54 


170565 


24 


79 7 


15 


52 13 


23 


64 12 


53 


163 4 


16 


44 13 


10 


32 41 


37 


111 4 


9 


30 2 


24 


75 


197417 


49 


181 3 


41 


114 6 


N/A 


N/A 





67 


228 1 


25 


96 9 


20 


65 19 


54 


185 2 


17 


59 4 


17 


52 


204367 


53 


181 4 


33 


86 6 


N/A 


N/A 





80 


249 1 


51 


143 5 


25 


72 16 


76 


236 3 


7 


18 2 


14 


46 


208217 


10 


34 7 


5 


17 12 


N/A 


N/A 





20 


69 4 


7 


29 13 


4 


15 37 


14 


70 6 


10 


33 6 


13 


40 



Table 5. Frequency analyses from cross correlations. Low-frequency peaks (0.4 mHz or below) were excluded from the noise estimates and not considered for 
significance. These 6 regions essentially cover: the spectrum bluer than the 6000 A gap, two regions redder than the gap and avoiding the weak telluric line 
region at AA 6275 — 6320 A, the Ha-region, the region with the Na D doublet, and a comparison region dominated by telluric lines. 
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a at 0.45 mHz, next-highest peak has A ma)[ = 7.7 ms 
b at 0.45 mHz, next-highest peak has A max = 8.6 m s~ 1 
c at0.48mHz, next-highest peak has A m n X = 31.3 ms - 1 



mine in a virtually approximation-free manner the mean magnetic 
field modulus (B) from measurement of the wave length separation 
of the resolved Zeeman components iMathvsll 19891) . (B) is the av- 
erage of the modulus of the magnetic vector, over the visible stellar 
hemisphere, weighted by the local line intensity. For a triplet pat- 
tern, its value (in G) is obtained from the wavelength separation 
between the central tt component and either of the a components, 
AA, by application of the formula: 

(B) = AA/(4.67 ■ 1(T 13 X 2 B g e g), (3) 

where A c is the central wavelength of the line and g c ff is the effec- 
tive Lande factor of the transition. Both AA and A c are expressed in 
A. For a doublet pattern, the relation between the field modulus and 
the wavelength separation of the split components (each of which 
is the superposition of a tt and a a component) is: 

(B) = AA/(9.34 • 1(T 13 X 2 C g cS ). (4) 

Note that Eq.[4] also describes the relation between (B) and the 
wavelength separation of the red and blue a components of a triplet. 

Mainly due to smearing by rotational broadening, only one 
of the studied stars (HD 208217) has clearly resolved lines from 
magnetic splitting that can be used with Eqs.[3]or|4]to estimate the 



field strength directly. In one other star, HD 107107, magnetic split- 
ting and rotational broadening are comparable for the lines with 
the highest magnetic sensitivity, so that it is possible to obtain an 
estimate of the surface magnetic field by fitting synthetic profiles 
to these lines. This estimate, which we shall denote by (B syn th), 
should be of the same order of magnitude as the mean magnetic 
field modulus, but it is not entirely equivalent to the latter. (In par- 
ticular, (B S y n th) is model-dependent, while (B) is not.) 

For the other stars we estimated the surface magnetic field 
from consideration of the magnetic broadening and intensification 
of magnetically sensitive lines. By assuming that a line's full-width 
at half maximum (FWHM) increases linearly with the separation of 
its (unresolved) Zeeman-split components, we used Eq.|4]to obtain 
an estimate (Bfwhm) of the magnetic field b y comparing lines 
with different effective Lande factors (see also iPrestonll 197 II) . The 
assumption is justified as long as the magnetic splitting of the anal- 
ysed lines is small compared to their overall width (in particular, 
due to rotational broadening); then, (Bfwhm) should typically be 
comparable to (B) or (B q ) (see below). 

Other estimates of the magnetic fields were obtained by com- 
paring the observed spectra to magnetically resolved or broadened 
line profiles to determine (B synt h) • These models were synthesised 
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with SYNTHMAG for a temperature grid of T cB = 7500, 8000 and 
8500 K. For this comparison, Cr and Fe lines with low Lande fac- 
tors were used to fix abundances and the rotation rate of the models 
prior to applying them to magnetically resolved or broadened lines. 
Instrumental broadening of 0.05 A was adopted while the macrotur- 
bulence was fixed for each star within the range 1 — 4 km s _ 1 . This 
approach worked well for magnetically broadened lines because 
both the equivalent width and the FWHM are affected by magnetic 
broadening; when absorption in a line is intrinsically spread over 
a larger wavelength range, the absorption becomes more effective, 
increasing the equivalent width (magnetic intensification). Again, 
the quantity that is derived should be comparable to the mean mag- 
netic field modulus, but in general is not exactly equal to it. 

Finally, the values of the mean quadratic magnetic field, 
(Bq) = ((B 2 ) + (B 2 )) 1/2 , were derived t hroug h application of 
the method described in iMathvs & Hubrid feOOfj) . Here (B 2 ) is 
the mean square magnetic field modulus (the average over the stel- 
lar disk of the square of the modulus of the field vector, weighted 
by the local emergent line intensity), while (B 2 ) is the mean square 
longitudinal field (the average over the stellar disk of the square of 
the line-of-sight component of the magnetic vector, weighted by the 
local emergent line intensity). The analysis was based on consider- 
ation of samples of reasonably unblended lines of Fe I and Fe II; 
the number of analysed lines varies from star to star and ranges 
from 14 to 33. This approach could not be successfully applied to 
HD 170565, for which the number of usable diagnostic lines (9) 
proved insufficient to untangle the contributions of the magnetic 
field (probably fairly strong), non-negligible rotation, and a possi- 
bly inhomogeneous distribution of iron on the stellar surface. The 
mean quadratic magnetic field is typically a few percent greater 
than the mean magnetic field modulus. 

For 3 stars, we also indicate in Table[2]values of the mean lon- 
gitudinal magnetic field (B z ) and the magnetic field modulus (B) 
from the literature. (_B Z ) is the average over the stellar disk of the 
component of the magnetic vector along the line of sight, weighted 
by the local emergent line intensity. This field moment strongly de- 
pends on the geometry of the observation, and contrary to, e.g., (B) 
or (Bq), it typically varies considerably during a stellar rotation cy- 
cle. Accordingly, it is not well suited to characterise the strength of 
the magnetic field in a star, other than to give (through its absolute 
value) a generally very conservative lower limit of the latter. But 
significant measurements of (B z ) at least provide definitive evi- 
dence that a star is magnetic. 

3.1.3 Radial-velocity shifts 

Precise radial velocity shifts were measured for several spectral 
lines with the centre-of-gravity method and by fitting with Gaussian 
profiles. We used the local continuum in the selected sub-regions of 
measured lines, which gave consistent amplitudes for both meth- 
ods in tests on UVES spectra for known roAp stars (see below). 
The centre-of-gravity method was preferred due to its stability and 
better ability to deal with blended line profiles. For strong and iso- 
lated lines, the two methods were comparable. Table[3] gives the 
most important spectral lines used, but the actual selection of lines 
depended on line blending and composition for each stellar case. 

Additionally, line shifts were determined by cross-correlating 
the spectra with averaged spectra of each series. The maxima of 
the correlation functions were determined with a spline fit, which 
was found more reliable and stable than with Gaussian or 4th- 
order polynomial functions. The cross-correlation regions were 
chosen to be free of static features: AA 5150 — 5800, 6035 — 6273, 
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Figure 3. Amplitude spectra for the known roAp star 33 Lib. Top: the 
2.015 mHz frequency for the Eu II 6645.06 A line. Bottom: Cross cor- 
relation (see text) for the region AA5150 — 5800 A. The frequencies 
2.015 mHz, its harmonic 4.030 mHz, and 1.769 mHz are recovered. Note 
different ordinate scales. 

6350 - 6700 and 6510 - 6570 A. However, to check for non-stellar 
periodicities, we also used two regions dominated by telluric or the 
interstellar lines of NaD: A A 5888 - 5898 and 6873 - 6899 A. 
Because the UVES pipeline produces merged spectra in the linear 
wavelength scale only, we rebinned the selected regi ons of the spec- 
tra to log wavelength which is more appropriate dTonrv & Davis! 
Il979l) when determining velocity shifts using large wavelength re- 
gions. For each region, the size of the bins was optimised according 
to pixel size and local spectral resolution. 

Because of the stratification of Ap star atmospheres, high- 
order mode pulsation amplitude and phase may change from ele- 
ment to element. Therefore the amplitudes of cross-correlation ve- 
locities cannot be directly compared to measurements of individual 
lines. Nevertheless, cross-correlation is efficient for detecting pul- 
sations by using long wavelength regions with numerous lines to 
obtain high S/N. 

IDL tools for line and cross-correlation measurements and 
analyses were developed and tested on UVES sp ectra of the known 
roA p stars 33 Lib and H P 154708, published bv lKurtz et al]fc005h 
and I Kurtz et alj d2006bl) . respectively. For 33 Lib, we confirm two 
frequencies, 2.015 mHz and 1.769 mHz (Fig.[3}, that are in ex- 
cellent agreement with the published frequencies, amplitudes and 
noise levels for lines of Eu II and Ha. In particular, we also con- 
firm the low-amplitude oscillation for Fe 1 5434.52 (31 ±4 m s _1 ). 
Similarly, we find amplitudes of 32 ± 5 ms _1 for Cal and Till 
lines and 130 ± 13 ms" 1 for Ball. A telluric line (6888.96 A) 
shows only noise with a highest amplitude of 9 ± 3 ms -1 . We 
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Figure 4. Same as Fig.[5] but for the known roAp star HD 154708. Top: 
significant amplitude is seen at the known frequency 2.088 mHz for the 
Prill 5299.99 A-line. Bottom: cross-correlation (AA5150 - 5800 A) un- 
ambiguously detects the 2.088 mHz frequency. 

also made cross-correlation measurements of the above-mentioned 
wavelength regions and recovered the 2.015 frequency with ampli- 
tudes in the range 20±2to50±5 ms" 1 . All cross-correlation and 
line measurements gave similar significance S/N = 10, although 
at very different amplitudes (see also Fig.[3}. 

The second test, the new roAp star HD 154708, was even 
stronger as this star pulsates with amplitudes tha t are among the 
smallest known for roAp stars. iKurtz et al.l d2006bl) needed to com- 
bine RV measurements for several lines of this star in order to de- 
tect its rapid oscillation. We confirm that no individual line (in- 
cluding Ha) shows signal on or above the 4 a level. Yet, as seen in 
Fig.|4] we directly recover the known 2.088 mHz mode with cross- 
correlations in the lower region of the spectra, AA5150 — 5800 A, 
with S/N — 6.2, and a marginal detection (110±32 ms _1 )using 
the single line Prill 5299.99 A. 

3.1.4 Frequency analysis 

Frequency analyses were performed u sing a Discrete Fourier 
Transform programm e dKurtzl Il985h and the PERIOD04 
dLenz & Bregeil 120051) programme. Linear trends in the indi- 
vidual ~ 2 hr radial- velocity series were fitted and removed with 
linear least-squares fitting. The noise a in the amplitude spectra 
(see Tables|4] and [5} is determined from lea s t-squ ares fitting 
of harmonics to the data following iDeemin j d 19751) . Because 
the barycentric velocity correction is approximately linear for 
each series of spectra, and rather small (the correction varies 
44— 150 ms -1 per hr for the 9 stars), it was eliminated with other 



drifts by a linear fit before the frequency analysis. With the 0.3 
arcsec slit and seeing conditions of 0.9 — 1.4 arcse c, the centring 
error for UVES is 50 — 100 ms -1 according to iBouchv et al.l 
d2004l) . Furthermore, a 1-mbar change in pressure may induce 
drifts of 90 ms -1 . During each of our observing nights, the 
pressure changed 1.5 — 2.0 mbar. We therefore expect instrumental 
drifts of up to 280 ms -1 per night, depending on seeing and 
pressure, and less during a 2-hr series on a star. The drift during 
a series of spectra may be non-linear (which is what we actually 
see in some cases). We noticed that comparison lines or regions 
with non-stellar constant lines occasionally exhibit drifts that are 
not seen for other regions of the same spectra (such as in Fig. |A10| 
panel 'Tell'). This difference may be because the strong and sharp 
telluric lines result in higher sensitivity, and they are influenced 
by fast wind speeds in the high layers of the Earth's atmosphere 
where telluric lines are formed. 

The null-results are presented on a star-by-star basis in Sec- 
tions l3.2H3.10l We emphasise the statistical fact that when calcu- 
lating about 50 amplitude spectra for each of nine candidate roAp 
stars, the chance for a spurious peak to reach the 4 a level in- 
creases. Furthermore, the analyses show several combined ampli- 
tude spectra with single, prominent peaks reaching 3 — 4<7. How- 
ever, the reality of such peaks is that they often originate from shal- 
low and blended lines. Our criteria for detecting rapid oscillations 
are therefore: i) a pea k of 4 times the noise (see lBreger et al J 19931 : 
iKuschnig et al.fl997r) in an amplitude spectrum, ii) confirmation in 
an amplitude spectrum for either a line, or combination of lines, 
of a different ionisation or element, or cross-correlation region, 
and Hi) only frequencies above 0.4 mHz (42min) are considered. 
The latter is because we have less control over drifts on these time 
scales in the wavelength calibration (i.e. we did not observe simul- 
taneous reference spectra). Lower frequencies are not typical for 
known roAp stars and may be due to, e.g., stellar rotation and sur- 
face spots. As an upper limit of the studied frequency range we use 
6 mHz. The sampling frequency is either 9.5 mHz or 15.6 mHz, and 
the Nyquist frequency is half of that but still above the frequencies 
in known roAp stars. 

In the following sections, we comment case-by-case on the 
individual stars in our sample. 

3.2 HD 107107 

With a magnitu de of V = 8.7 34, this star is one of the faintest 
in our sample. iMartinezl d 19931) obtained 1.89 h photometry dur- 
ing a single night and excluded periodic variability above 0.3 and 
0.7 mmag for frequencies higher than 1 .0 and 0.4 mHz respectively. 
Based on his f3 ph otometry (Tabled the corr esponding temperature 
from the grids bv lMoon & Dworetskvl dl985l) is 8300 K. The Hip- 
parcos mission obtained 100 useful measurements; an amplitude 
spectrum of those has a noise level of 3.1 mmag and excludes peaks 
above 9 mmag. The distribution of Hipparcos data is, however, not 
well suited for detecting rapid oscillations, but may instead show 
rotation periods of spot ted stars. The Michigan Spectral Catalogue 
(Houk &Cowlevll 19751) classifies the star as Ap CrEuSr. No spec- 
troscopy has been previously published for HD 107107. 

Our observations of this star comprise 111 UVES spectra ob- 
tained over a time span of 1.94 hr with 63-s time resolution. The 
spectra are rotationally broadened to vsini = 10.5 kms" 1 (Ta- 
ble[2jl and have single and strong REE lines. The star is very pe- 
culiar with strong Nd III and Pr III lines, probably with a large ratio 
between abundances of singly- and doubly-ionised REEs. The lat- 
ter could indicate ionisation disequilibria of REEs, a common fea- 
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Figure 5. Amplitude spectra from cross-correlation for the HD 107107 
spectra. Top: for the wavelength region AA5150 - 5800 A; Bottom: for 
the wavelength region AA 6350 - 6700 A. 



ture among known roAp stars (Rvabchikova et al. 2004i). The NaD 
AA 5889.95 and 5895.92 A lines each have a stellar and two inter- 
stellar components that are sharper, stronger and red-shifted with 
respect to the stellar one. 

Radial velocity shifts of 47 stellar lines were measured (Ta- 
ble[4]> using the full line profiles where line blending permitted it. 

The core of Ha is constant to 65 m s" 1 (a = 20 m s" 1 ) and other pulsationally stable above 39 m s~ (a = 12 m s" ) for all Nd 



Figure 6. Magnetic broadening measurements of 25 Fe lines of HD 107107. 
Measured Gaussian FWHM (AA) vs product of laboratory wavelength 
squared and Lande factor (A^ g c g ). A significant relation (r = 0.77) fitted 
with a least-squares linear fit, is indicated with a line. 



the cross-correlation radial velocity shifts is around 4 ms' 1 with- 
out significant peaks above 15 ms _ . This is comparable to the 
result for the stable telluric line region (AA6873 — 6899 A). 

Some Cr II and Fe II lines are magnetically broadened or par- 
tially resolved and the mean quadratic field determined for 27 
iron lines is (73 q ) = 5.2 ± 0.4 kG. Using measured FWHM 
of 25 Fe lines, we derive (Fig.[6jl a field strength of 3.4 ± 
0.6 kG (1 a error). Due to the combination of marginally resolved 
lines, rotational broadening and line blending, the use of ap- 
parently double lines gave less consistent results that, however, 
support existence of a strong magnetic field. SYNTHMAG syn- 
thetic profiles for different magnetic field strengths were computed 
and compared to Cr II 5 116.049, Cm 53 18.382, Cr 1 5247.566 and 
Fe II 6149.25 A. The best fit was obtained for a magnetic field 
strength of (-B sy nth) = 4.5 ± 0.5 kG (estimated error). The mag- 
netic field modulus was further estimated directly to be (B) = 
5.6 ± 2.3 kG from separations of the partially resolved components 
of 5 Cr I, Cr II and Fe I double lines using Gaussian fitting. 

HD 107107 is, therefore, a chemically peculiar A star that is 



lines also exhibit no detectable variability in the considered fre- 
quency range (0.4 — 6.0 mHz). Tableland Figure |A3l show selected 
results of the atomic lines measured in the frequency analysis. As 
indicated in the table, radial velocity s eries for different l ines of 
same species were combined following iKurtz etal (2006a) to re- 
duce the noise. Combining all line measurements, including dif- 
ferent species, we reach a — 8 ms" 1 and a maximum amplitude 
of 34 ms" 1 . One peak at 1.8mHz (panel 'HD 107107 all') is just 
above the 4 a detection limit, but originates from weak and blended 
Ndlll lines. Yet, this ion has some of the highest amplitudes in 
roAp stars and is excellent for detecting rapid pulsations. However, 
this peak cannot be confirmed by lines of other elements, so follow- 
ing the criteria in Sect. l3.1.4l we rather label it a possible detection 
that needs confirmation. Other lines (such as Eu II) show peaks near 
the detection limit, but again the frequencies are unconfirmed else- 
where. This demonstrates the difficulty in reliable detection of sig- 
nal at this noise level. Selected results of the cross-correlation anal- 
ysis, presented in Tableland Fig.|5] strengthen this null-result. For 
the 5 stellar wavelength regions defined in Sect. l3.1.3l the noise in 



and Pr lines combined, and 9 ms" 1 (a = 3 ms" 1 ) for cross- 
correlations. The star is a new magnetic star, with marginally re- 
solved Zeeman-split lines and a field of (B q ) = 5.2 ± 0.4 kG. 



3.3 HD 110072 

This is the faintest star in our sample. iMartineJ i 19931) excluded 
photometric variability above about 1.4 and 0.8 mmag for frequen- 
cies above 0.4 and 0.9 mHz respectively, based on 67 min of pho- 
tometry on a single night. The Hipparcos data show no s ignifi- 
cant p eaks above 22 mmag at shorter frequencies. iHouk & Cowlevl 
d 19751) correctly note that HD 1 10072 is type Ap Sr(Cr) rather than 
K0 dSAO Staffll966l classification source: M. W. Mayall). 



We collected 69 spectra with UVES in 2.06 hr with a time res- 
olution of 107 s. The spectra are sharp lined: the estimated v sin i 
is only 3.3 ± 0.5 kms" 1 and many strong REE lines are visible, 
e.g., Pr, Nd, Y, Eu II and Ce II. Also Cr, Fe II, Ni II, Co I and Al II 
are strong. In addition to a stellar component, the Na D lines each 
have 3 sharper and bluer (with respect to the stellar component) in- 
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terstellar components of comparable strengths. The average spec- 
trum is very similar to those of two known roAp stars 33 Lib and 
HD 176232 (Fig.|A]}. These have temp eratures of T eS = 7550 ± 
150 and 7550 ± 100 K respectively jRvabchikova et all l2004h . 
which is supported by the indistinguishable shapes of the Ha wings 
of 33 Lib and HD1 10072 ( T cff = 7300 K). All three stars have sim- 
ilar peculiarities for REEs, and Ba II, Si I and Ca I are considerably 
weaker in HD 110072. 

The Ha profile is strongly asymmetric (see Fig.[2} with a dip 
~ 70 kms -1 blueward of the Ha core. This dip is only 3 — 4 per 
cent below the Ha wing, but about 40 per cent broader (FWHM) 
than the Ha core. We re-observed the star 2 yr later with FEROS 
at the ESO 2.2-m telescope and found this feature to have dis- 
appeared. This indicates that HD 110072 is a binary with a sec- 
ondary star that may be less luminous than the primary. The broad- 
ness of the core of the secondary's Ha line can either be due to 
faster rotation (~ 50 kms -1 ) or a later spectral type (Ha weak- 
ens toward GO). At other wavelengths, the UVES and FEROS 
spectra are largely identical and no secondary spectrum is seen 
(see also Fig. lAlb . However, a few lines appear only in the re- 
cent FEROS spectrum (resolution R = 48000) such as at locations 
of Sc 1 615 1.20 A, 1 6156.77 A, Smil 6157.53 A, Nd II 6549.52 A 
and Fe II 6552.33 A. These lines have a broadening similar to the 
primary's spectrum. We suspect they originate from this star and 
their appearance is a result of viewing different aspects of a mag- 
netic field and/or spotted chemical surface distribution combined 
with line blending. 

Radial velocity shifts were measured for 49 stellar lines (Ta- 
bleO and show no detectable signal. The blended Ha core is stable 
to 135 ms~ but the noise of a = 52 ms" 1 in this radial veloc- 
ity series is considerable, caused by fewer spectra and lower S/N 
than obtained for the other roAp candidates. Highest significance 
(3.9 cr) is seen for the combined Fell lines at 1.86 mHz, but iron 
is known in other roAp stars to have low amplitudes, and there 
is no support for this frequency from other lines. Combining all 
49 lines ("Fig. lA4b reduces the noise to 8ms -1 , excluding peaks 
above 25 ms~ , The cross-correlations in Table[5] with examples 
in Fig.[7] result in flat amplitude spectra above 0.45 mHz. A 4 a 
peak at 0.45 mHz is caused by non-linear drifts. The telluric line 
region is stable to 17 ms -1 . 

The mean quadratic field determined for 33 iron lines is 
{Bq) = 1.5 ± 0.6 kG. Comparison of 14 Cr and Fe lines with 
SYNTHMAG models indicates broadening by a magnetic field of 
1 — 3kG, in particular for lines such as Cr II 5246, Crl5247, 
Fe 1 5324 and Fe II 6149.25 A. Measurements of FWHM of 27 iron 
lines with IRAF's onedspec.splot task give (Fig.[8jl a mean magnetic 
field modulus of (Bfwhm) = L5 ± 0.3kG (1 a error). 

HD 110072, the coolest star in our sample, with T e g = 
7300 K, is thus a sharp-lined, double-lined binary with strongly 
peculiar lines of, e.g., Nd, Pr and Eu. HD 110072 is a new mag- 
netic star with a field of (B q ) = 1.5 ± 0.6 kG and has spec- 
tral features very similar to the roAp stars 33 Lib and HD 176232. 
It is therefore intriguing that the star is pulsationally stable to 
42 ms -1 (a = 10 ms -1 ) for all Nd and Pr lines combined, and 
12 ms -1 (a = 3 ms -1 ) for cross-correlations. The slow rotation 
of HD 1 10072, its peculiar abundances and the rare combination of 
a magnetic field and its binary status makes it an important case for 
studying stellar evolution and diffusion processes. 
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Figure 7. Amplitude spectra from cross-correlation for the HD 110072 
spectra. Top: for the wavelength region AA5150 - 5800 A; Bottom: for 
the wavelength region AA 6350 - 6700 A. 
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Figure 8. Magnetic broadening measurements of 27 Fe lines in HD 1 10072 
plotting the measured Gaussian FWHM (AA) vs the product of laboratory 
wavelength squared and Lande factor (Aj? g e g)- A significant relation (r = 
0.70), fitted with a least-squares linear fit, is indicated with a line. 
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Figure 9. Amplitude spectra from cross-correlation for the HD 131750 
spectra. Top: for the wavelength region AA5150 - 5800 A; Bottom: for 
the wavelength region AA 6350 - 6700 A. 



3.4 HD 131750 

For this star, Houk & Cowlevl fl97^ 's classification is Ap CrEuSr. 
Martin^ i 19931) observed it 1 — 2h on each of three nights. The 
first night showed a clear peak around 8.5 mHz (which is outside 
the range we consider), but the other nights showed no variabil- 
ity abo ve 0.6mmag for / > 0.6 mHz. Istrohmeier. Fischer. & Ottl 
( 1966) listed an unconfirmed 0.35 mag variability from photo- 
graphic plates (no period given). Hipparcos data do not support 
this. 

Our 111 spectra, obtained in 1.98 hr with a time resolution 
64 s, show HD 131750 to be rotating with v sin i = 25.3 kms" 1 
with strong lines of Nd, Eu and Pr. The NaD doublet is strong 
with multiple sharp interstellar components. Radial velocity shifts 
of 40 stellar lines were measured (Table|4j- The noise in the radial- 
velocity measurements is relatively high due to the rotation rate. 
The core of Ha is stable to 72 ms" 1 (a = 21 ms" 1 ), while all 
measured lines combined (Fig. |A5t reveal no rapid pulsation above 
39 ms -1 (a = 8 ms" 1 ). The cross-correlations provide an upper 
limit of 9 — 15 ms -1 to the averaged radial velocity shifts of the 
measured regions. Flat amplitude spectra of two of these regions are 
given in Fig.|9] The NaD region has a significant peak at 0.45 mHz 
caused by non-linear instrumental drifts in the data. Elsewhere, the 
amplitude spectra are flat, devoid of significant peaks. 

Some of the stellar line profiles have flat 'squared' cores, such 
as Cm 5613.18 and Fe II 5854.19, and others have apparently split 



cores such as Fell 5457.73, Cm 5472.60 and Fe 1 5862.35. Mag- 
netically resolved lines require a Zeeman splitting at least compa- 
rable to rotational broadening, in this case around (B) = 8 kG. 
The mean quadratic field was derived to (B q ) — 5.3 ± 3.3 kG by 
using 9 iron lines. It is our impression that the upper limit on any 
magnetic field modulus (at this rotation phase) is 8 kG. 

HD 131750 is thus an Ap star with strong REEs, no rapid os- 
cillations above 58 ms" 1 (a — 18 ms -1 ) amplitude for all Nd 
and Pr lines combined, and 9 ms" 1 (a = 3 ms" 1 ) for cross- 
correlations. The rotation rate v sin i = 25.3 km s" 1 results in con- 
siderable line blending and the star may have a magnetic field of a 
few kG. 

3.5 HD 132322 

iHouk & Cowlevl dl975h classify this star as 'Ap SrCrEu, AT and 
note that Sr is extremely stron g. Photometric Stromgren indices 
also indicate strong peculiarity. Martinet i ll 9931) examined the star 
during a single night and excluded rapid photometri c variability 
above 0.4 mmag for the frequency range 0.4— 10 mHz. lLevato et al.l 
4l996h obtained 2 spectra of the star and found v sin i = 85 kms , 
a som ewhat high value for an Ap star. Additionally, iHubrig et al.l 
J2006h discovered a mean longitudinal magnetic field of 357 ± 
51 G. 

HD 132322, V = 7.357 mag, is the second brightest star in 
our sample and the individual spectra have S/N well above 100. 
A total of 1 1 1 spectra were obtained in 2.03 hr, providing a 66-s 
time resolution. The spectra appear to show splitting of all lines, 
such as Ba, Fe and Cr lines, and also of all REEs, such as Nd and 
EuII. Surprisingly, even lines with small Lande factors are dou- 
ble (Fig. 1 10b which would indicate a spotted surface distribution 
rather than splitting of lines by a magnetic field. It is, however, im- 
probable that all elements, in particular Fe, are split due to a spot- 
ted surface. In the course of the 2 hr of observations, we do not 
find any systematic change in the radial velocity difference of the 
components, nor in their centre of gravity. The estimated rotational 
broadening is v sin i = 25-35 km s~ 1 , based on the full double pro- 
files, and v sin i = 34.3 ±1.8 km s" 1 from the quadratic magnetic 
field analysis, which together with the peculiarity results in consid- 
erable blending. Lines of Pr are not very strong but may also be 
double. The otherwise useful lines of Cell are absent. For some 
roAp stars this ion has the highest pulsation amplitude, such as for 
/3CrB (Kurtz, Elkin & Mathys 2007). The NaD doublet is strong 
with a broad and a sharp (interstellar) component. 

Radial velocity shifts measured for 34 stellar lines put an up- 
per limit of 24 ms" 1 (a — 8 m s" 1 ) to rapid pulsation when com- 
bining all lines (Fig. lA6b . The core of Ha is stable to 59 ms -1 
(<r = 19 ms -1 , Table|4]l. All amplitude spectra are flat, except 
for that of combined Cm which has a non-significant peak at 
1.31 mHz. The cross-correlation analysis (see, e.g., Fig.ll It shows 
flat amplitude spectra down to amplitudes of 6 — 15 m s" 1 without 
any significant peaks. 

Using 14 iron lines, the mean quadratic field could only be 
constrained to a maximum intensity of (B q ) ^ 6.0 kG. This weak 
constraint results from the fact that, because of the considerable 
broadening and distortion of the spectral lines, only a small num- 
ber of them could be identified as sufficiently free from blends to be 
used in the analysis. Indeed many of the lines, including those that 
are magnetically insensitive, show double structures, which typi- 
cally consist of a component separated 30 ± 3 km s~ 1 from a red- 
der and 45 ± 15 per cent weaker component (in equivalent width). 
The FWHM of the weaker component is 31 ± 18 per cent less than 
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Figure 10. Double line structures in the magnetically insensitive line 
Fe 1 5434.52 in the averaged spectrum of HD 132322 (thick line). A syn- 
thetic model for vsini = 30 kms" 1 and (B) = 0G is superposed (thin 
line) with its dominant atomic lines indicated. 



for the other. This pattern is similar for lines of REEs, Cr and Fe (cf . 
Figs.|A2]and Fig.[l0]l. A comparison of the Ha profile to a synthetic 
spectrum, shifted in wavelength corresponding to the separation in 
the double structures, firmly excludes a secondary spectrum of a 
star of comparable brightness as it would have introduced a strong 
asymmetry. New high-resolution spectra are required to understand 
these double structures. 

HD 132322 is a magnetic Ap star with projected rotation ve- 
locity v sin i = 34 km s" 1 . All lines, even those of REEs, are dou- 
ble except for Ha. This can partly, but not fully, be explained by 
abundance spots or a secondary spectrum. Pulsations are excluded 
down to 37 ms" 1 (a — 13 ms" 1 ) for all Nd and Pr lines com- 
bined, and 6 ms -1 (a — 2 ms -1 ) for cross-correlations. 



3.6 HD 151301 



The star is cla ssified Ap SrCrEu by iHouk & Cowlevl Jl975h . 
Martinez (1993) observed it on 6 nights for 0.9 — 1.4 hr each, and 
excluded pulsations down to 0.5 — 0.8 mmag for frequencies above 
0.5 mHz. 

We obtained 111 UVES spectra in 1.95 hr at a time resolution 
of 63 s. The lines of REEs are strong (Nd, Pr, Eu) and many are 
double (see also Fig.[2} indicating abundance spots. The NaD lines 
have a stellar component and a sharp, strong, interstellar line at 
longer wavelength. The photometric temperature is T e g — 8000 K 
and the spectroscopy provides an upper limit T e « = 9000 K for 
log g = 3.5. With the astrometric luminosity of the star, this agrees 
with HD 151301 being more than halfway through its main se- 
quence lifetime. 

Radial velocity shifts of 39 stellar lines were measured (Ta- 
ble[4]l. Lines of PrII, Cell and Ndll are weak with consider- 
able scatter in their radial-velocity series. When combining all 
39 lines, we find an upper limit on rapid pulsation of 18 ms" 
(a — 6 ms" 1 ), while the core of Ha is stable to 54 ms -1 
(a = 19 m s _1 ). Each double component of the two strongest Eu II 
lines was measured and the combined amplitude spectrum (Panel 
'EuIF, Fig. |A7| l shows a significant (4.2 cr) peak at 2.22 mHz. It is 
equally significant for the individual Eu II lines and when combin- 
ing three available Eu II lines. No other line or combination of lines 
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Figure 11. Amplitude spectra from cross-correlation for the HD 132322 
spectra. Top: for the wavelength region AA 5150 — 5800 A; Bottom: for the 
wavelength region AA 6350 - 6700 A. 



confirm this peak, including a few weak lines of La II and Ce II 
added to the analysis. Even cross-correlation of the whole EuII 
profiles did not recover the 2.22 mHz frequency and it is therefore 
considered a probable, but unconfirmed detection. The spectrum 
of HD 151301 is rich in lines, and the cross-correlations reach a 
low noise level (Tabled, in particular in the bluer spectrum below 
the 6000A gap (1.4 ms" 1 , Fig.l 12t. The spectrum region above 
the 6000 A gap results in considerably larger noise (8-10 ms" 1 ) 
but also excludes significant rapid pulsations. The rotational broad- 
ening is vsini= 13.7 ± 1.1 kms" 1 . Some Fel and Fell lines 
show additional broadening or asymmetric profiles. An upper limit 
of the mean quadratic field is found to (B q ) ^ 2.4 kG, using 
18 Fe lines. SYNTHMAG fitting to Cr 115313.56, Fei6137.69 and 
Fel 5266.55 rejects magnetic fields stronger than 2 kG. However, 
from line- width measurements of magnetic broadening of 13 Fe 
lines, a weak, possibly insignificant relation (r = 0.50) indicates 
(Bfwhm) = 1.2 ± 0.60 kG. 

HD 151301 isastrongly chemically peculiar star possibly hav- 
ing a magnetic field of up to 2 kG. The surface distribution of REEs 
is spotted. No pulsations are seen down to amplitudes of 24 ms" 1 
(a = 8 ms -1 ) for all Nd and Pr lines combined, and 5.4 ms" 1 
(a — 1.4 ms -1 ) for cross-correlations. 



3.7 HD 170565 

Martinez (1993) observed this star on 5 nights. In general, no pul- 
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Figure 12. Amplitude spectra from cross-correlation for the HD 151301 
spectra. Top: for the wavelength region AA 5150 - 5800 A; Bottom: for the 
wavelength region AA 6350 — 6700 A. 

sations are seen down to 0.7 mmag for frequencies above 0.5 mHz. 
However, the last two nights show two peaks around 1.7 and 
2.6 mHz at 3 — 4 times the noise. The first of th ese is also present 
on the first and most intensively observed night. iKudrvavtsev et al.l 
d2006h detected a magnetic field in this star with a mean longitudi- 
nal field of 1.76 ± 0.17 kG. 

We obtained 85 spectra of the star in 2.50 h with a time resolu- 
tion of 106 s. The star is highly peculiar and has v sin i = 18 ms -1 . 
There are many strong REE lines, such as those of Eu II, Ce II and 
Nd that are all double, indicating abundance spots. PrII is less 
strong, but also double. As a result line-blending is considerable. 
Radial velocity shifts of 41 stellar lines were measured and re- 
ject pulsations down to 32 ms -1 (a = 8 ms" 1 ) when combin- 
ing all lines (Fig. lA8l and Table|4]l. The core of Ha shows stabil- 
ity to 75 ms" 1 (a — 24 ms -1 ). There are no confirmed signifi- 
cant peaks in periodograms for any combination of radial velocity 
series. Cross-correlations show stability down to 10 — 30 ms -1 
(Tabled and Fig.QJJ. Only the AA6350 - 6700 A region shows 
a significant (4.2 cr) peak at 2.59 mHz (the same as the second 
'transient' photometric period). No other line measurements or 
cross-correlations support this unconfirmed detection. Cr is dou- 
ble or broadened while Fe is broadened in several cases. Magnetic 
measurements were tried with three methods: mean quadratic field 
measurements, SYNTHMAG fitting to 4 Cr and Fe lines, and from 
widths of 22 Fe lines. However, due to the rotation, line blending, 
and a possible inhomogeneous stellar surface distribution of iron, 
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Figure 13. Amplitude spectra from cross-correlation for the HD 170565 
spectra. Top: for the wavelength region AA 5150 - 5800 A; Bottom: for the 
wavelength region AA 6350 - 6700 A. 



it is not possible to constrain the known magnetic field with these 
data. 

This star is a known magnetic star which is consistent with the 
present data. It is a chemically peculiar Ap star with a spotted sur- 
face distribution of Cr and REEs, and is pulsationally stable down 
to 32 m s~ 1 (a — 10 m s~ 1 ) for all Nd and Pr lines combined, and 
10 ms -1 (a — 5 ms -1 ) for cross-correlations. 



3.8 HD 197417 

The c lassification for this object is Ap CrEu(Sr) faouk & Cowlevl 
and the star has been investigated earlier due to its chemical 
peculia rities and photom etric variability associated with its rotation 
period. Marti ne?fl993h observed the star on two nights. The cor- 
responding amplitude spectra are fla t and place an upper limit of 
0.6 mmag on photometric variability. iFloquet et alj d 19841) studied 
this star and found a photometric rotation period of 4.551±0.002 d. 
With spectra obtained at 12 A mm - dispersion, they used Balmer 
and Call lines to determine T c g = 9500 K and \ogg — 4.0, 
and from Mg II 448 1 A they measured v sin i = 23 km s~ 1 . A strong 
variability was noticed in the intensity and profile of Ca II 3922.6 A 
and by assuming an oblique rotator geometry and a spotted sur- 
face distribution of elements, they proposed an inclination of i — 
64° (angle of the stellar rotation axis to the line of sight) and that 
Ca II, Eu II and Sr II in particular seemed located in a common spot. 
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Figure 14. Amplitude spectra from cross-correlation for the HD 197417 
spectra. Top: for the wavelength region AA 5150 - 5800 A; Bottom: for the 
wavelength region AA 6350 — 6700 A. 

From 156 Hipparcos measurements we do not find this period sig- 
nificant, but cannot reject i t based on the noisi er Hipparcos data 
alone. Based on 4 spectra, iLevato et alj Jl996h listed the star as 
proba bly single (25 per ce nt chance for random velocity distribu- 
tion). [PMnzen~eTaO |2005|) list the star as a confirmed chemically 
peculiar star with Aa = 0.054 and (b - y) Q = 0.032. 

HD 197417 has the second-highest Sci (0.015) of our sample, 
and the 125 UVES spectra, obtained in 2.20 h at 63-s time resolu- 
tion, show weak lines of REEs; PrII, Prill and, e.g., Nd III 6550.32, 
6145.07 and 6327.24 A are almost absent (2 — 4 per cent below the 
continuum). Both Eu II 6437.68 and 6645.06 A lines are weak (less 
than 5 per cent below the continuum). Some REE lines are partially 
split, e.g. Prill 6866.80, Euli 6437.64 and Eu 11 6645.06 A, which 
may indicate that the stellar surface is spotted. This may in com- 
bination with the rotation v sin i = 25.5 ms -1 partly explain the 
weak REE lines of this Ap star. Lines of Co I, Fe II and La II are 
strong, while a weak Ca I supports the high photometric tempera- 
ture, T cB = 8400 K. Ce II is absent. The Na D lines have a stellar 
component and a sharper interstellar component at longer wave- 
lengths. 

Radial velocity shifts of 33 stellar lines were measured (Ta- 
ble|4] and Fig. |A9t . The noise is considerable due to the low pe- 
culiarities and line blending due to rotational broadening, but 
the radial-velocity series for the core of Ha alone, or all 33 
lines combined, exclude rapid pulsation to 48 ms -1 (a — 13). 
Cross-correlations result in flat amplitude spectra down to 5 — 



10 ms 1 (see, e.g., Fig.l 14b. Magnetically sensitive lines, such as 
Cr 115116.04, Fe 116149.25, and Fel6232.64A, are not magneti- 
cally resolved. Comparison of 8 Fe and Cr lines with magnetically 
broadened SYNTHMAG profiles, and also an attempt to measure the 
mean quadratic field using 20 lines, exclude magnetic fields above 
~ 2 kG. Due to considerable line blending, smaller fields cannot be 
quantified without polarimetric measurements. 

HD 197417 does not appear to have strong REE abundances in 
the wavelength regions we cover. This may, however, be an effect 
from a spotted surface distribution of REEs and rotational broad- 
ening. The existence of an inhomogeneous surface distribution of 
REEs is supported by the star's known (rotational) photometric 
period and the many double REE lines. A spotted surface on Ap 
stars can often be related to presence of a magnetic field, for which 
we in this case put an upper limit at 2kG. The star is stable to 
65 ms -1 (a = 20 ms -1 ) for all Nd and Pr lines combined, and 
5 ms - (cr = 2 ms -1 ) for cross-correlations. The significant ro- 
tation v sin i = 25.5 kms -1 does, combined with blending from the 
many double lines, limit the lines available for analysis. 

3.9 HD 204367 

The Michigan Spectral Catalogue lists the star as A(p SrEuCr) and 
notes that it is either a weak Ap star, or a normal star of spectral type 
AOIV/Vs. Supposedly due to this, iMartinezI d 19931) did not make 
any time series photometry of this object. The star has the highest 
Sci (0.020 mag) and Co index in our sample. Because HD 204367 
is one of the least evolved of the studied stars (Fig.[T|l, the high 
indices rat her indicate a less pecu liar spectrum rather than a higher 
luminosity. Man froid et alj dl998) obtained 21 measurements (one 
per night) in the Geneva photometric system of this Ap candidate 
and found no variability (mmag level) over a 23-night run. With 90 
Hipparcos measurements, we find the star stable to 6.5 mmag for 
periods longer than a day. 

We obtained 111 spectra of this moderate rotator 
(usini= 10.8 kms -1 ) in 1.96hr at a time resolution of 64s. 
The spectra revealed the least chemically peculiar star in this study, 
with considerably fewer lines and longer continuum windows than 
any of the other stars. Lines of Nd, Pr, Ce and Cr are either absent 
or much weaker than in the other studied stars. The absorption of 
Eu II is in all cases less than 2 per cent below the continuum and 
only a few REE lines could be used in the velocity analysis. Lines 
of Ni, Si, S, Na, Ba and Zn are strong. 

Radial velocity shifts of 38 stellar lines were measured. The 
noise is considerable, but when combining all lines (Fig. lAlOl l, we 
can exclude rapid pulsation with amplitudes above 36 ms -1 (a — 
12 ms -1 ), and 46 ms -1 (a = 14 ms -1 ) for the Ha line core 
alone (Table|4]l. There are no significant peaks in the periodograms. 
Cross-correlations produce similar 'flat' amplitude spectra (Table|5] 
and Fig.ll5t. An upper limit of the mean quadratic field is found 
with 27 Fe lines to <B q ) < 1.4 kG, while SYNTHMAG models 
compared to 16 Cr and Fe lines gave an upper limit of (B syn th) ^ 
1 kG. Further, line width measurements of 33 iron lines showed a 
weak relation (Fig,|16t that would indicate a rather weak field of 
0.77 ± 0.24 kG (1 a error). This fit is, however, barely significant 
as a significance test only gave t = 2.9 while t > 3.0 is required. 
Polarimetry is therefore needed for verification. 

This star shows no detectable periodic variability down to 
72 ms" 1 (o — 25 ms -1 ) for all Nd and Pr lines combined, and 
5 ms -1 (a — 2 ms -1 ) for cross-correlations. Rotation is slow 
(v sin i = 10.8 km s~ ), the spectrum is nearly devoid of REEs, and 
the star may have a small magnetic field of ~0.8 kG. 
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Figure 15. Amplitude spectra from cross-correlation for the HD 204367 
spectra. Top: for the wavelength region AA 5150 - 5800 A; Bottom: for the 
wavelength region AA 6350 — 6700 A. 
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Figure 16. Magnetic broadening measurements of 33 Fe lines of 
HD 204367. Measured Gaussian FWHM (AA) vs product of laboratory 
wavelength squared and Lande factor (A;? g a g). A weak relation (r = 0.50) 
fitted with a least-squares linear fit, is indicated with a line. 
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Figure 17. Amplitude spectra from cross-correlation for the HD 208217 
spectra. Top: for the wavelength region AA 5150 — 5800 A; Bottom: for the 
wavelength region AA 6350 - 6700 A. 



3.10 HD 208217 

Thi s known magnetic and peculiar star, classified as Ap SrEuCr, 
Al dHouk & Cowlevl[l 975) has been previously examined in sev- 
eral studies. iMartinej dl993l) spent nearly 13 h on the star on 8 sep- 
arate nights. With Scj — —0.19 it is strongly peculiar. The high- 
speed photometry excludes periodicities above 0.4 mmag (frequen- 
cies above 0.8 mHz) and 0.8 mm ag (frequencie s abov e 0.5 mHz). 
With spectroscopic observations, iMathvs et al.l (l997) detected a 
mean magnetic field modulus varying with a semi-amplitude of 
nearly 1000 G about a mean value of 7958 G (a = 588 G) over a 
rotation period of 8.44475 d dManfroid et aljl997h . For this period, 
the epoch of our spectroscopic observations occurs near a negative 
extremum of the longitudinal field (Mathys, unpublished observa- 
tions). For roAp stars, maximum pulsation amplitude occurs when 
one of the stellar magnetic poles come into sight. In this case the 
negative magnetic pole, and our chances for detecting roAp pulsa- 
tions should therefore be optimal. Mathys et al. furthermore found 
the star to be a single-lined binary with a most likely orbital period 
of at least 2 yr. 

With a typical S/N of 150 per spectrum, this set of 138 UVES 
spectra obtained in 2.52 h with a 66-s time resolution has the high- 
est quality of this study. The spectra show strong lines of Nd, Pr 
and Eu II, while Ca I is weak which could indicate a higher temper- 
ature than the T e g — 7700 K from photometry, but the wings of Ha 
agree well with models for T e ff in the range 7500 — 8000 K. Lines 
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Figure 18. Magnetic field strength of HD 208217 demonstrated with the 
average UVES spectrum (thick line) for the absorption lines Fe I 5434.52 
(top panel) and Fell 6 149.25 A (bottom) which respectively have a low 
and a high Lande factor. Two SYNTHMAG models are superposed for 
mean field moduli of 0.0 (full line) and 8.0 kG (dashed line). The mod- 
els use i)sini = 9.4kms _1 , [Fe/N] = -5.4 (Fel5434.52) and -4.4 
(Fell 6149.25). 



of REEs such as Nd II 6549.52 are mostly single; Eu II 6437.64 A is 
double while the single Eu II 6437.64 appears broadened. The lines 
of Ce II are absent. 

Radial velocity shifts of 52 lines were measured (Tabled and 
show for all lines combined no rapid oscillations above 14 ms" 1 
(cr = 4 ms -1 ). The core of Ha is stable to 40 ms -1 (a = 
13 ms -1 ). However, lines of Prill and EuII show peaks of 3.9 a 
and 4.9 a respectively but for different frequencies and are re- 
jected. The cross-correlations produce flat amplitude spectra reach- 
ing the lowest noise level in this study: the integrated radial- 
velocity noise is only 1.4 — 2.8 ms - and no peaks appear above 
3 a. Lines of REEs are in general single, however Eu II 6437.64 is 
broadened while EuII 6437.64 is double. Also Fell 6149.25 and 
Fell 6432.68 A are double, and iron lines are generally broad- 
ened, consistently with the strong magnetic field known for this 
star. The rotation is low (t>sini = 10.8 kms -1 ), but the spectrum 
is distorted by the magnetic field and blends from the abundant 
peculiar elements, so line identification and analysis is compli- 



cated. Using 16 Fe lines, the mean quadratic field is found to 
(Bq) = 8.0 ± 0.5 kG. With SYNTHMAG fitting of 17 lines of dif- 
ferent Lande factors (see examples in Fig.l 1 8t the field strength is 
found to be (B S ynth) = 7.5 ± 0.5 kG, while the Zeeman splitting 
of Fell 6149.25 and Cr II 5318 A gives a magnetic field modulus of 
8.0±0.3kG. 

Summarising, this star has strongly peculiar REE abundances, 
a (known) strong magnetic field of {B q } = 8.0 ± 0.5 kG and is sta- 
ble to 15 ms -1 (a — 4 ms -1 ) for all Nd and Pr lines combined, 
and 3 ms -1 (cr = 1 ms -1 ) for cross-correlations. 



4 DISCUSSION 

The class of roAp stars is notoriously difficult to supplement with 
new members, as demonstrated by several photometric and spec- 
troscopic studies before this one. But motivate d by the recent dis - 
covery of the luminous roAp star HD 116114 dElkin et af]|2005l) . 
we have spectroscopically tested a sample of 9 luminous Ap stars 
for rapid pulsations. Using lines known to show pulsations in roAp 
stars, we reach typical upper amplitude limits in radial velocity of: 
40 - 75 ms" 1 (cr = 13 - 24 ms" 1 ) for the line core of Ha, 
20 — 65 ms" 1 (a — 7 — 20 ms" 1 ) when combining all measured 
Nd and Pr lines, and 20 — 40 ms -1 (cr = 7 — 11 ms -1 ) when 
combining all measured lines. With cross-correlations, using large 
wavelength regions, we typically reach upper amplitude limits of 
4—10 m s -1 (a = 1 — 4 ms" 1 ). In spite of a clear theoretical 
prediction dCunhall2002h and empirical (HD 116114) evidence for 
roAp pulsations in this part of the Hertzsprung-Russell diagram, 
we end up with 9 null-results, or noAp stars. A number of ques- 
tions are therefore pertinent to discuss. 

How well does our test sample resemble known roAp stars ? 

All studied stars have strong REE lines, except for HD 204367 and 
possibly also HD 197 '417 '. They have the core-wing anomaly typical 
for roAp stars, and Hipparcos luminosities with our temperature es- 
timates place them inside the predicted roAp instability strip. Sev- 
eral appear to have spotted surface distributions of REEs (such as 
HD 170565, HD 151301, HD 132322 and perhaps also HD 197417) 
which is typically associated with the strong magnetic fields com- 
mon in known roAp stars. Indeed most of the stars are magnetic, 
and cover a range in magnetic field strengths of 0.4 — 8 .0 kG, com- 
parable to that of known roAp stars dKurtz et al. 2006b). In the case 
of the sharp-lined HD 110072, we compared its spectrum in detail 
with those of two known roAp stars, and found remarkable sim- 
ilarities for the REEs. However, HD 110072 and HD 208217 are 
double-lined and single-lined binaries, respectively, which might 
indirectly influence their stability to high-f requency pu lsations by 
reducing the magnetic field intensity (see [Cunha 2002] and refer- 
ences therein). Still, the orbits are probably too wide in both cases 
for tidal interaction to occur and HD 208217 has a known strong 
magnetic field ((B) = 8 kG). The known roAp stars have cases 
of wide binaries, such as /3CrB (spectroscopic binary), HR3831, 
aCir, 7Equ and HD 99563 (visual binaries). In these regards, the 
studied sample has the characteristics of roAp stars. 

Could pulsations have been overlooked? 

Kurtz et al. (2006b) published radial velocity amplitudes for the 
Ha cores of 16 roAp stars. Of these, only 3 have amplitudes be- 
low 75 ms" 1 (HD 116114, HD 154708 and HD 166473, of which 
two have amplitudes above 3 cr), while the rest range from 148 — 
2528 ms -1 . Further, radial velocity series for Pr and Nd line 
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measurements in UVES spectra of these 16 roAp stars (Kurtz 
et al., partly unpublished), show that about 75 per cent of the 
measured and significant (3 a) amplitudes are within the range 
350 — 1600 ms _1 . Five of these stars have very small Nd and 
Pr amplitudes (60 - 90 ms" 1 ): HD 166473, HD 116114, /3CrB, 
33 Lib and HD 154708. In such difficult cases, other lines or combi- 
nations of several lines makes detection of pulsations possible. We 
successfully tested our procedures on the two latter roAp stars in 
Sect. |3.1.3l and also used combinations of several lines, including 
of different elements. More of the other 19 known roAp stars have 
low amplitudes, such as lOAql, but our tests and analyses show 
that we reach these amplitude levels and should have detected such 
rapid pulsations if present in the studied sample. 

A complication for our analyses is typical v sin i ~ 10 — 
30 km s _1 combined with double lines due to either spots on the 
stellar surface and/or magnetic splitting, which results in consider- 
able line blending and makes line identification and analysis more 
difficult. However, our radial-velocity analysis was based partly on 
cross-correlations that are more robust than line measurements (as 
shown by our tests for two roAp stars) and this method similarly re- 
sults in flat amplitude spectra that exclude rapid oscillations to rel- 
atively small roAp- amplitude levels. It also seems improbable that, 
e.g., unfavourable viewing angles of the global pulsations or short 
mode lifetimes could explain a momentary lapse of detectable pul- 
sation amplitudes simultaneously in all nine stars. In fact, we know 
independently that HD 208217 was observed near its magnetic neg- 
ative extremum where roAp pulsations are expected to have maxi- 
mum amplitudes. Future surveys like this one may benefit from be- 
ing repeated at different rotation phases. We also note that the near- 
normal REE abundances of HD 197417 and HD 204367 reduce the 
probability that they are roAp stars, given the strong peculiarity of 
all known roAp stars. 

Do these stars really not pulsate ? 

In pulsators such as roAp stars oscillations are intrinsically unsta- 
ble. Their excitation depends on the balance between the driving 
and damping of the oscillations over each pulsation cycle. In roAp 
stars this balance is thought to be particularly delicate. On one 
hand, the amount of energy input through the opacity mechanism 
acting on the hydrogen ionisation region depends strongly on the 
interaction between the magnetic field and envelope convection, 
being m aximal in the regions where envelo pe convection is sup- 
pressed l lBalmforth etalTl200ll: ICunhj|2002h . On the other hand, 
the direct effect of the magnetic field on pulsations can introduce 
significant energy losses, through slow Alfven waves in the interior 
and through acoustic waves in the atmosphere, bo th resulting from 
mode conversion in the magnetic boundary layer dCunha & Goughl 
2000; Saio 2005). Due to this delicate balance, it is not too sur- 
prising that roAp and noAp stars occupy the same locus in the HR 
diagram. Despite the developments in theoretical studies of linear 
non-adiabatic pulsations in models of roAp stars, we still lack a 
theoretical study that takes into account all these phenomena si- 
multaneously and, thus, cannot firmly predict the conditions under 
which puls ations should be expec t ed in roAp stars. In fact, both 
studies of lCunhj j2002h and Saio d2005h considered the extreme 
case in which envelope convection is fully suppressed. Moreover, 
the first of these studies did not consider the direct effect of the 
magnetic field on pulsations and neglected the energy losses as a 
result of mode conversion, and the second study, while consider- 
ing mode conversion, assumed the waves are fully reflected at the 
surface, hence neglecting energy losses through acoustic running 
waves in the atmosphere. 



As discussed by ICunhal d2002l) . the condition for suppres- 
sion of envelope convection, which seems necessary to make the 
high frequency modes unstable, is in principle harder to fulfil in 
evolved stars due to the increase with age of the absolute value of 
the buoyancy frequency in the region where hydrogen is ionised. 
Hence, it is likely that in evolved stars oscillations are excited 
only if the magnetic field is relatively strong. Unfortunately, the 
complexity of the interaction between magnetic field and convec- 
tion makes it impossible to derive a global convective stability 
criterion, even if local criteria fo r convective stability m ay be es- 
tablished dGough & Tavle3ll966l: iMoss & TavleJll9691) (see also 
iTheado et alj 120051 , for an extensive discussion on this subject). 
Thus, the magnetic intensity needed to suppress convection at a 
given age, for a given mass, is very hard to establish. The more 
evolved roAp star HD 1 161 14, in which a relatively low frequency 
oscillation was found well in agreement with theoretical predic- 
tions, has a magnetic field modulus of ~ 6 kG. In contrast with 
this, most stars in our sample have estimated mean magnetic field 
moduli around or below 2 kG. The clear exceptions are HD 107 107, 
HD 131750 and HD 208217. Of these three, the latter is clearly an 
important test case to check the theoretical predictions. It has the 
strongest confirmed magnetic field in our sample and is strongly 
peculiar. However, we observed HD 208217 when one of its mag- 
netic poles were almost visible, so pulsation should have been near 
its maximum amplitude. 

From the observational point of view, one way to investigate 
the conditions under which roAp star oscillations are excited, and 
thus test theoretical models, is by identifying systematic differences 
between roAp and noAp stars. This study would have been able 
to detect pulsations in all the known roAp stars, and any missed 
rapid pulsations must have amplitudes lower than these. Hence we 
conclude that based on the obtained data, most stars in our sample 
are indeed noAp stars, and also excellent roAp candidates. Despite 
this conclusion, it is premature to state that we can confirm the 
evidence that noAp stars are in average more luminous and more 
evolved than roAp stars, as indicated by earlier studi es based on 
photometric surveys for pulsations in roAp candidate s dNorth et"al] 
1 19971 ; lHandler & Paunzerj|l999l ; iHubrig et al.ll2000h . To conclude 
that, we would also need to search for rapid pulsations in a con- 
trol group of less evolved roAp stars with lower luminosity and 
compare the frequency of null-results in the two cases. Such a sur- 
vey has been started and results are expected in the near future. 
The next step is to analyse the noAp stellar atmospheres in detail, 
taking temperature gradients and the abundance stratification into 
account. 

New spectra are needed of HD 132322 to clarify the origin of 
its double line structures, and of HD 1 10072 to verify its secondary 
spectrum and to test for the 'spurious' absorption lines in the recent 
FEROS spectrum. Polarimetry of HD 110072 and HD 204367 is 
needed to confirm the detected magnetic fields. 
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Table Al. Abundances in log N/N to t for our model spectrum used for 
line identification. Effective temperature T e g = 8750 K and surface 
gravity is log g = 4.0. 
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A2 Selected regions of average spectra 




Figure Al. Comparison of the averaged spectrum of HD 1 10072 (full line) with averaged UVES spectra of the two known roAp stars HD 176232 (dot-dashed) 
and 33 Lib (offset and dotted). Note the similarities in strong REEs (Nd III and Pr III), while Ba II, Si I and Ca I all are considerably weaker in HD 1 10072. The 
spectrum of the fainter HD 1 10072 was smoothed over every 4 pixels. Note the absence of lines from a secondary at this wavelength region of HD 1 10072. 
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Figure A2. Comparison of the profiles of all examined Ap stars. To increase the readability, the individual spectra have been offset in flux, and appear with 
increasing temperature upwards. The selected region shows examples of double lines due to magnetic splitting, such as HD 208217, or spots on the surface of 
the stars. The spectrum of HD 1 10072 was smoothed over every two pixels. 
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A3 Amplitude spectra for combined line measurements 

These figures show star-by-star amplitude spectra calculated for 
combined radial- velocity series for elements as indicated by top 
labels: 'all' indicates all available lines, 'NdPr' all lines for 
neodymium and praseodymium. For centre-of-gravity line 
measurements. Note that for double lines, both components were 
measured separately and included. Compare with Table|4] 

This paper has been typeset from a TjSC/ ETgX file prepared by the 
author. 



24 L. M. Freyhammer et al. 



HD107107 Window 



HD107107 Fel 



HD107107 Fel 



m 150 - 
F 



2 3 

f mHz 



HD107107 Ndl 









\A rA lA- 







2 3 

f mHz 



HD107107 Prl 




HH 1 07 1 07 Prill 



b 150 - 
E 



2 3 

f mHz 



HD107107 NdPr 



i u,..o. 

2 3 4 

f mHz 

HD107107 Cell 






f mHz 
HD1 071 07 Crll 



2 3 

f mHz 



HD107107 oil 






1 2 



3 

f mHz 



HD107107 Nd 



g 100 
E 



. j..vji... y., v . .yy.w.vy^. 



2 3 

f mHz 



HD107107 Pr 



.5 100 
E 




' » Wo/ V 



2 3 

f mHz 



HD107107 Eul 



g 100 

E 




Figure A3. Amplitude spectra for HD 107107. Panel labels indicate the element lines(s) in the combined RV series. Top left panel is the window function for 
a 3 mHz signal. Panels 'all' and 'Tell' are, respectively, for all available lines combined and for all telluric lines combined. 



A search for rapid pulsations among 9 luminous Ap stars 25 




Figure A4. Same as Fig. lA3l but for HD 1 10072. The Nyquist frequency is 4.6 mHz. 
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Figure A5. Same as Fig.|A3]but for HD 131750. 
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Figure A6. Same as Fig. lA3l but for HD 132322. The Nyquist frequency is 4.7 mHz. 
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Figure A7. Same as Fig.|A3]but for HD 151301. 
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Figure A8. Same as Fig. lA3l but for HD 170565. The Nyquist frequency is 4.7 mHz. 
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Figure A9. Same as Fig.|A3]but for HD 197417. 
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Figure A10. Same as Fig.|A3]but for HD 204367. 
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Figure All. Same as Fig.|A3]but for HD 208217. 



